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Abstract 

We investigate different neutrino signals from the decay of dark matter particles 
to determine the prospects for their detection, and more specifically if any spectral 
signature can be disentangled from the background in present and future neutrino 
observatories. If detected, such a signal could bring an independent confirmation 
of the dark matter interpretation of the dramatic rise in the positron fraction 
above lOGeV recently observed by the PAMELA satellite experiment and offer 
the possibility of distinguishing between astrophysical sources and dark matter 
decay or annihilation. In combination with other signals, it may also be possible 
to distinguish among different dark matter decay channels. 



1 Introduction 



While the existence of dark matter is now firmly established [T], little is known about 
the properties of the particles that make up the dark matter, including their lifetime. 
The dark matter particles are often assumed to be perfectly stable as the result of 
a symmetry, e.g. i?-parity in supersymmetric models. However, from the gravitational 
evidence for the existence of dark matter we can only infer directly that the dark matter 
has to be stable on timescales comparable to the age of the Universe. Measurements of 
cosmic-ray antimatter, gamma rays and neutrinos, on the other hand, typically impose 
much more stringent constraints on the dark matter lifetime in the respective decay 
channels [2]. 

As with baryonic matter itself, there are good reasons to consider the case of unstable 
dark matter. From the theoretical point of view, we expect at least gravity to violate 
any global symmetry (and in some cases the symmetry breaking takes place at a lower 
scale, as e.g. in models of Grand Unification |3]). Therefore, we can expect the presence 
of at least non-renormalisable operators in the theory allowing for dark matter decay. In 
other cases, the dark matter parity symmetry may be only approximate from the start or 
spontaneously broken, as it occurs in some models of -R-parity breaking [4] . Alternatively, 
the coupling involved in the decay may be very strongly suppressed, as in the case of 
a tiny kinetic mixing between visible sector and hidden sector [5]. In those cases it is 
natural to expect a very long lifetime for the dark matter particle, which may exceed the 
age of the Universe by many orders of magnitude. Nevertheless, even for such extremely 
long lifetimes the decay signals may be in the observable range pi IH El El [71 El El CDl E] • 

Another important clue in relation to decaying dark matter is the observation of 
several excesses in the fluxes of cosmic rays in the energy region above a few GeV. Namely, 
the PAMELA experiment observed a steep rise in the positron fraction extending up to 
at least 100 GeV [12]. Furthermore, the experiments Fermi LAT, H.E.S.S. and ATIC 
measured the total e'^e~ flux, finding that this flux is harder than expected and falls 
off steeply above 1 TeV [I3]. Such signals can be well explained by the decay of dark 
matter particles in the Galactic halo if the decays are sufficiently "leptophilic," such that 
one may avoid the overproduction of antiprotons. In fact, the observed antiproton flux 
is in complete agreement with the flux expected from production by spallation of cosmic 
rays [H] . Regarding the mass of the decaying particle, the PAMELA data require it to be 
larger than 200 GeV since the rise in the positron fraction extends up to at least 100 GeV, 
while the Fermi and H.E.S.S. measurements of the total e~^e~ flux indicate a mass of a 
few TeV, depending on the background [15] and the particular decay modes. The electron 
flux from dark matter decay required to fit the data points towards a lifetime of around 
lO^^s in all cases [7]. 
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Similar signatures can also be generated by annihilating dark matter, but in this case 
some tension with the constraints coming from the measurement of radio emissions from 
the centre of the Galaxy [16] and from inverse Compton scattering with starlight and 
the cosmic microwave background [IT] is present. Note that astrophysical explanations 
for the positron excess have also been put forward [18], the most popular being one or 
several nearby pulsars as additional sources of electrons and positrons. 

The decay of dark matter particles is an interesting and viable explanation of the 
observed electron anomalies, and it is worthwhile trying to confirm or exclude this possi- 
bility by a complementary examination of other indirect detection channels, in particular 
neutrinos [SlEj. Neutrinos, from this perspective, have two clear advantages. Firstly, they 
are unaffected by magnetic fields and thus, like photons, allow to reconstruct the direc- 
tion of their origin; therefore they would offer a clear way to distinguish between the 
cases of annihilating and decaying dark matter, as well as pulsar interpretations of the 
signal (even assuming that the pulsars produce also neutrinos in the energy range consid- 
ered). Secondly, they are typically produced along with or from the decay of the charged 
leptons in many "leptophilic" decaying dark matter models [10]. In such cases, therefore, 
the flux of neutrinos is correlated with the other cosmic-ray signals, and their spectrum 
may give direct information on the dark matter decay channel. In particular, choosing 
the mass and lifetime of the dark matter particle such as to yield a good agreement with 
the PAMELA positron excess, one can directly predict the rates for the corresponding 
neutrino signal and look for it in present and future experiments. These are the two 
advantages we will try to exploit in this paper. 

On the other hand, neutrinos also suffer from some clear disadvantages with respect 
to other indirect detection channels: The large atmospheric neutrino background makes 
it difficult to disentangle any signal up to TeV energies for the lifetimes indicated by the 
cosmic-ray anomalies mentioned above. A further disadvantage is the necessity of very 
large detectors to measure the comparably small neutrino fluxes expected from dark 
matter decay. Fortunately, new large neutrino detectors, namely IceCube and possibly 
KMSNeT, will become fully operational in the near future and may allow to detect even 
the small signals we discuss here. 

This paper is organised as follows: In section [2] we will discuss the neutrino flux 
expected from decaying dark matter in our Galaxy and compare it with the one from 
dark matter annihilation in order to discuss the best strategy for the detection of the 
signal in these two cases. In section [3], we will present the spectral signatures for a number 
of different dark matter decay modes. In section H] we will give the present bounds from 
neutrino experiments and the expected rates for present and future neutrino detectors. 
We will also discuss the prospects for distinguishing between different neutrino spectra 
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in case a signal is detected. We will finally present our conclusions in section [51 



2 Neutrino Fluxes 



We concentrate in this section on the neutrino flux expected from the dark matter in the 
Milky Way halo, since on one hand it is the dominant source, and on the other hand it 
has a nontrivial directionality that may be exploited, as in the case of gamma rays [19], to 
disentangle the different hypotheses of dark matter decay versus annihilation. In addition, 
an isotropic extragalactic component is expected from unresolved cosmological sources, 
which in the case of decaying dark matter is of the same order of magnitude as the 
halo contribution, so that it may increase the overall signal by a factor of two or so. 
This extragalactic component is expected to be negligible in the case of dark matter 
annihilation. 



dE 47rrDM"ZDM dE 

l.o.s. 



Decaying Dark Matter For the case of decaying dark matter particles in the halo, 
the differential flux of neutrinos is given by the following integral along the line of sight: 

phMljdl, (2.1, 

where tdm and ttidm are the lifetime and the mass of the decaying particle, is the 
neutrino energy spectrum from the decay and Phaio is the dark matter density in the halo. 
Adopting for phaio the NFW density profile, we obtain for the averaged full-sky flux 

d_J^ \ . pU^o+iJ ^ r ( ^\ _ _ 1 1^2 ^ 

dE I 47rrDMmDM dE \ V + ^J 2 \ 

= 1.3 X 10"* (cm^ s sr)~^ — ^ , (2.2) 

V ^DM / V "^DM / dE 

with the local halo density pioc = 0.3GeVcm~^, the solar distance from the Galactic 

centre i?© = 8.5 kpc and = 20kpc. The numerical result is only weakly dependent on 

the halo parameters and the profile. 

The flux is inversely proportional to the product of the dark matter particle mass and 
lifetime. Thus, for a fixed lifetime the flux is inversely proportional to the dark matter 
mass due to the lower number density of dark matter particles for higher masses. 



Annihilating Dark Matter For an annihilating particle, the differential flux of neu- 
trinos is instead given by the following integral along the line of sight: 

^■^halo _ {(^v)dm dK f 2 /A .r 

~dE~ ~ ~dE J ' ^^-^^ 



l.o.s. 



where {o'v)j^j^ and ttidm are the dark matter annihilation cross-section and the dark 
matter mass, is the neutrino spectrum from annihilation instead of decay, and the 
line-of-sight integral, in this case, contains the square of the halo density. It is there- 
fore clear that for the same halo profile the annihilating dark matter signal is strongly 
enhanced towards the centre of the Galaxy, especially for cuspy halo profiles. 

Propagation After the neutrinos are produced in the decay or annihilation of dark 
matter particles, they travel in straight lines through the Galaxy, essentially without 
any interactions. The only modifications to the fluxes during this time are due to flavour 
oscillations [20]. In fact, using the experimental best-fit values for the neutrino mixing 
angles, sin^ 6'i2 = 0.304, sin^^23 = 0.5 and sin^ 6*13 = 0.01 [21], and neglecting possible 
CP-violating effects, the neutrino oscillation probabilities in vacuum are given by 



Thus, a primary neutrino flux in a specific flavour is redistributed almost equally into all 
neutrino flavours during propagation and any flavour information is lost. On the other 
hand, this means that nearly the same signal is present in any flavour and may allow to 
choose the best channel for discovery according to the background and efficiency of the 
detector. 

2.1 Background Fluxes 

Let us now discuss the background for our neutrino signal. The main background for the 
observation of neutrinos in the GeV to TeV range are neutrinos produced in cosmic-ray 
interactions with the Earth's atmosphere. Here we use the atmospheric neutrino fluxes 
calculated by Honda et al. [22]. The theoretical uncertainty of these fluxes is estimated 
to be better than 25 % in the GeV to TeV range, while the uncertainty in the ratio of 
the different flavours is significantly smaller. We extend the atmospheric neutrino fluxes 
to energies above 10 TeV using the slopes given by Volkova et al. [23] . 

Conventional electron and muon neutrinos are directly produced from pion and kaon 
decays. While electron neutrinos are practically unaffected by neutrino oscillations due to 
the large oscillation length, muon neutrinos, particularly at low energies, can be converted 
into tau neutrinos and provide the dominant tau neutrino background at energies below 
1 TeV. The conversion probability of muon neutrinos into tau neutrinos is given by 



P(z/e ^ Z/e) 



0.56 



0.39. 



0.22 



(2.4) 




(2.5) 
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In this expression, Eiy is the neutrino energy and L is their propagation length after being 
produced in the atmosphere, which is given by 



with ~ 6.4 X 10^ km being the Earth's radius and /i ~ 15 km the mean altitude at 
which atmospheric neutrinos are produced. 

In addition to the conventional atmospheric neutrino flux from pion and kaon decays 
there is a prompt neutrino flux from the decay of charmed particles that are also pro- 
duced in cosmic-ray collisions with the atmosphere. The prompt neutrinos have a harder 
spectrum than the conventional ones and therefore dominate at higher energies (roughly 
above 10 TeV for electron neutrinos and above 100 TeV for muon neutrinos). Since these 
contributions are not well understood and in any case subdominant in the energy range 
that is of interest here, we neglect them in the present study. 

On the other hand, the prompt tau neutrinos start to dominate around 1 TeV (and 
at even smaller energies for downgoing neutrinos). Thus we include this contribution 
using the parametrisation [23] 



where x = log^o (S(GeV)), A = 6.69, B = 1.05, C = 0.150 and D = -0.00820. This 
parametrisation is valid in the energy range of 100 GeV up to 1 PeV. However, we point 
out that compared to the conventional atmospheric neutrino flux the prompt flux suffers 
from larger uncertainties. 

Other neutrino backgrounds in the considered energy range are neutrinos produced 
in cosmic-ray interactions with the solar corona [25j and those produced in cosmic-ray 
interactions with the interstellar medium in the Milky Way [26]. While the former is 
subdominant in diffuse searches for all flavours [8] and can be excluded from the analysis 
by excluding neutrinos from the direction of the Sun, the latter represents an irreducible, 
ill-understood neutrino background for searches in the Galactic disc direction. In fact, the 
flux of Galactic neutrinos is expected to become comparable to the atmospheric electron 
neutrino background for the Galactic disc direction and energies in the TeV range. 

2.2 General Detection Strategy and Use of Directionality 

In view of the subdominant neutrino signals from dark matter decays it is important to 
devise strategies that reduce the background. In [S] it was proposed to use directionality 
in order to reduce the background in the tau neutrino channel. This is possible since the 
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Figure 1: Significance of the signal as a function of tlie cone half angle towards the Galac- 
tic centre normalised to the significance of the full-sky signal for annihilating/decaying 
dark matter depending on the different density profiles. 

tau neutrino background at low energies comes mainly from the muon neutrino oscillation 
and is therefore strongly suppressed in the zenith direction. Some of the authors of [9] 
propose instead to search for an enhanced muon neutrino signal only in the direction of 
the Galactic centre. However, taking into account the typically low neutrino event rates, 
it is not always the best strategy to optimise the signal-to-background ratio. Instead, the 
statistical significance a = S/\/B (number of signal events divided by the square root of 
the number of background events) is a better measure for comparing different detection 
strategies. 

In Figure [1] we show the significance of the signal as a function of the cone half angle 
around the Galactic centre normalised to the significance of the full-sky observation. 
Here we assume a background that (practically) does not depend on Galactic coordinates 
like atmospheric neutrinos and neglect the contribution of Galactic neutrinos, which is 
subdominant in the case of muon neutrinos. We see clearly that for annihilating dark 
matter the best way to detect the signal is indeed looking towards the Galactic centre: 
The cone half angle offering the best signal-to-square root of background ratio varies 
depending on the cuspiness of the profile, but it is always between ~ 0° (NFW) and 30° 
(isothermal). Note that in any case the gain of looking at the Galactic centre is not very 
large for a cored profile like the isothermal one. 
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For the case of decaying dark matter on the other hand, the best strategy is to 
measure the full-sky signal and not concentrate on the region around the Galactic centre. 
In fact the gain coming from the enhanced dark matter density is counteracted by the 
smallness of the collecting area and so the significance of the signal goes quickly to zero 
as a function of the angle for any profile, even for cuspy profiles like the NFW profile. 
The observation of only a fraction of the sky around the Galactic centre direction leads 
to an increase in the signal-to-background ratio, but not of the significance. We therefore 
conclude that for decaying dark matter there is no advantage in looking only at the 
Galactic centre. The full-sky signal offers not only better statistics, but also a higher 
significance. 

Considering the directionality of the atmospheric background instead of the signal, 
another good strategy might be to exploit the fact that the flux from the zenith direction 
is (dependent on the energy) a few times smaller than from the horizontal direction. 
Assuming a signal that does not depend on the zenith angle, the observation of only a 
fraction of the sky around the zenith direction is again clearly leading to an increase in 
the signal-to-background ratio. Also in this case though, it turns out that the best value 
for the significance is achieved for a full-sky observation. 

We can therefore conclude that exploiting the directionality of the signal or back- 
ground, apart from the case of specific flavours like the tau neutrino discussed in |8], 
is not promising for the first detection of decaying dark matter. The largest rate and 
significance is achieved for a full-sky search, and this is the option we will discuss in the 
following. On the other hand, directionality offers a clear way to disentangle decaying 
dark matter from either annihilating dark matter, where looking into the Galactic centre 
should give an increase in significance, or from point sources like dwarf galaxies, pulsars 
and supernova remnants. 

3 Neutrino and Muon Spectra 

The neutrino spectra depend on the decay channel of the dark matter particle. The 
simplest possibility is a direct decay into two neutrinos for a scalar particle or into 
Z^v for a fermion. Then the resulting spectrum is just a monochromatic line for the 
Galactic signal and an integral of the redshifted line from the extragalactic signal (and 
a continuum contribution from the fragmentation of the boson in the case of decay 
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into Z'^u). So for this case we have the simple spectra 



V 



(DM oc 5 E 



^DM 

2 



) 



dE 



+ Ci.o.s. \ \^^^ 
(+ continuum) , 




where X = z/ and we have assumed here that the mass of X is neghgible. Ci.o.s. 

is the ratio of the extragalactic and the Galactic signal, which is a number of order one 
given by 



and K, = Qj>^/Qrn ~ 3 is the ratio of the dark energy density and the matter density in 
the Universe. 

Another characteristic spectrum is that obtained from a three-body decay into three 
leptons which has the familiar triangular shape when plotted on a logarithmic 

axislll In this case the expression for the extragalactic signal is more involved, but it still 
appears as a softer triangular shape which dominates at low energies, as can be seen 
from the change in slope for the three-body spectra in Figures [2] and [31 We show here as 
an example the spectrum of the decay of a fermionic dark matter particle mediated by 
a heavy scalar particle, corresponding to a scalar-type 4-fermion interaction. 

Finally, a continuous neutrino spectrum is generated by any heavy particle that 
decays into neutrinos, like the muon, or fragments into charged pions, like the electroweak 
gauge bosons or the tau lepton. We use here as examples of continuum neutrino spectra, 
the spectra arising from the decay of a scalar particle into longitudinal gauge bosons 
similar to the Higgs decay and from the decay of a fermionic chiral lepton into W^i^, Z^u. 
We use PYTHIA 6.4 [27] to simulate the gauge boson fragmentation and the heavy 
leptonic decays. 

The corresponding neutrino fluxes for these different types of spectra are shown in 
Figures [2] and [3] for a scalar and a fermionic dark matter candidate, respectively, together 
with the expected atmospheric background and the data measured by the Frejus [28] . 
Super-Kamiokande [29], AMANDA-II [30] and IceCube [3l] experiments. We see that 
for a lifetime of the order of 10^^ s, which is the order of magnitude suggested by the 
PAMELA excess [11], the signal always lies below the measured background of muon 

^This holds for the most common scenarios, where the decay is mediated by a heavy scalar or a heavy 
vector boson. In both cases the Michel parameter p is equal to 3/4, yielding the same neutrino energy 





spectrum. 
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Figure 2: Neutrino spectra for different decay channels of a scalar dark matter candidate 
compared to the expected background of atmospheric neutrinos from Honda et al. [22] 
and the data of Frejus [28], Super-Kamiokande [29] . Amanda- 1 1 [30] and IceCube [31]. The 
flux is computed for a dark matter mass of 1 TeV {top) or 10 TeV {bottom) and a lifetime 
of 10^^ s. The line from the two-body decay into uu and the extragalactic contribution to 
this decay spectrum is easy to distinguish. The spectra from the decays of a dark matter 
candidate into t~^t~ , Z^Z^ or W^W^ are softer at the endpoint. The low-energy 

tail of these decay channels is due to the muon/tau decay and Z^ jW^ fragmentation. 
Due to the steeply falling atmospheric background the signal-to-background ratio at the 
endpoint of the decay spectra increases significantly for larger dark matter masses. 



10 



atmospheric neutrinos 



DM ^Zv 
DM ^ eev/nnv/xTV 
DM ^ We/Wn/WT 
Super-K v„ 
Amanda-ll v 
Frejus Vg 
Frejus v„ 
Amanda-ll v 
lceCube-22 v„ 




10" 



10"' 



w 

> 

CD 

S 10-5 

> 

LU 

P -fi 

^ 10^ 

T3 

X 

LU 10 



atmospheric neutrinos 



DM ^Zv 
Dl\|/I -> eev/^nv/xTV 
D VI ^ We/Wn/Wx 
Super-K v„ 
Amanda-ll v 
Frejus Vg 
Frejus v„ 
Amanda-ll v 
lceCube-22 v„ 



10"' 




10-" 
Ev (GeV) 



Figure 3: Neutrino spectra for different decay channels of a fermionic dark matter can- 
didate compared to the expected background of atmospheric neutrinos from Honda et 
al. [22] and the data of Frejus [28], Super-Kamiokande [29], Amanda-II [30] and Ice- 
Cube [31]. The flux is computed for a dark matter mass of 1 TeV {top) or lOTeV {hot- 
torn) and a lifetime of 10^^ s. The pure line and three-body decays are easy to distinguish 
and correspond to the line contained in Z^u and the decay into e^e~v. Note that the 
low-energy tail of Z^v and of the other leptonic three-body decays is due to the Z*^ frag- 
mentation and muon/tau decay. Also shown are the cases of a pure continuum spectrum 
coming from the decay into W^i^. Due to the steeply falling atmospheric background 
the signal-to-background ratio at the endpoint of the decay spectra increases significantly 
for larger dark matter masses. 
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neutrinos. The best signal-to-background ratio is achieved for the high-energy end of the 
spectrum, which gives information about the mass scale of the decaying particle. 

The neutrino spectra shown in Figures [2] and |3] look rather distinctive, and an inter- 
esting question is whether they can be disentangled in a neutrino detector. We have to 
consider that neutrino detectors do not really measure neutrinos directly, but the corre- 
sponding charged leptons or showers, produced in the interactions of neutrinos with the 
intervening matter. 



3.1 Neutrino Interactions 

Since we are interested in neutrino energies much larger than nucleon masses we only 
take into account deep inelastic neutrino-nucleon scattering. Neutrino-electron elastic 
scattering is subdominant in this energy range and will be neglected. 

The cross-sections for deep inelastic scattering of (anti)neutrinos off nucleons at rest 
are given by 



dy TT 



\ '^CC/NC + "CC/NC \^ ~ y) 



3.2 X 10"^« ^ (a-;^, + 6-;^, (1 - yf 



(3.3) 



with flee" ~ 0.15, 0.25, 6^0" ~ 0.04, 0.06 and a^^'" = &cc^' ^cc" ~ ^cc^ charged- 
current interactions, and a'C,^"' = 0.058, 0.064, fe^c" ~ 0.022, 0.019 and aj^^" = 



ivp,n up,n 



for neutral-current interactions [201 [32] . The inelasticity y is given by 

y = l-— or y^ -— , (3.4) 

where is the energy of the generated lepton and i^^had is the energy of the generated 
hadronic shower. Equation (13.31) holds only for neutrino energies up to the TeV region 
when the effect of the massive gauge boson propagators cannot be neglected anymore. 
For higher energies the cross-sections are overestimated. 

For the total neutrino-nucleon cross-sections one obtains 

"cc/Ncy^^> — ^ -^f I "cc/NC ^ 3 "cc/NC ) • y^^-"^' 

As we can see, the total cross-section is proportional to the energy of the incoming 
neutrino in the considered energy range. 



3.2 Muon Neutrinos 

The charged- current deep inelastic scattering of a muon neutrino off a nucleus pro- 
duces a hadronic shower and a muon. These track-like events can be clearly identified in 
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material density (g/cm^) {Z/A) a (GeVcm^/g) /3 (cm^/g) 



water 



standard rock 



ice 



2.650 



0.918 



1.000 



0.55509 



0.55509 



0.5 



2.3 X 10"^ 
2.7 X 10-3 
2.7 X 10"3 



4.4 X 10"' 
3.3 X 10-' 
3.3 X 10-' 



Table 1: Density, proton- number-to- mass-number ratio and approximate muon energy 
loss parameters for the materials of interest in Cherenkov detectors. 

Cherenkov detectors via the Cherenkov light cone of the relativistic muon. 
3.2.1 Through-going Muons 

Since muons are rather long-lived (cr^ = 658.650 m), their range is only limited by energy 
loss during their passage through matter and not by their lifetime. Therefore Cherenkov 
detectors can also observe muons that are generated in the surrounding material of the 
detector. This effect enhances the effective detector area for high-energy muon neutrinos. 

The average rate of muon energy loss can be written as 



where a{E^) describes the ionisation energy loss and P{E^) takes into account the energy 
loss due to radiative processes: e^e~ pair production, bremsstrahlung and photonuclear 
contributions. Both a^E^) and /3(-E'^) are slowly varying functions of the muon energy. 
As long as we can approximate a and /3 as energy- independent, the average range after 
which the muon energy drops below a threshold energy Ejj^ is given by 



where p is the density of the medium. The relevant parameters for standard rock, water 
and ice are given in Table [TJ The values of the density and the average proton-number- 
to- mass-number ratio are taken from [33] . The muon energy loss parameters given in the 
table are best-fit values from the fit of equation ( 13. 7p to the tabulated data in [33] . From 
equation ( 13. 7p we can determine the initial muon energy as a function of the final muon 
energy and the muon range: 



In fact, equation ( 13. 6 p does not account for the stochastic nature of radiative muon energy 
losses which start to dominate at TeV energies {E > a//3), and therefore equation (13. 7p 
overestimates the muon range for large energies. 




(3.6) 




(3.7) 
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The rate of muon neutrino induced through-going muon events is given by 

oo 



dN 



dn / dE„ 



dQ / dE„ 



dE, 



dJuf,{Eu^,(), 



dE, 



da^^c{E,^,E,) 



dE„ 



np + {p n) 



(3.9) 



M 



xR,{E,,El')A^;{E,,e, 



where the number density of protons is given by rip = pN^ (Z/A) and the density of 
neutrons by n,„ = pNa{1 — (Z/A)). = 6.022 x lO^^mol"-'^ is the Avogadro constant, 
p is the density of the material and (Z/A) is the average ratio of the proton number 
and the mass number of the material as given in Table [H Due to the small neutrino- 
nucleon cross-section the attenuation term that accounts for the absorption of part of 
the signal and background neutrino fluxes during the passage of the Earth is negligible 
in the considered energy range. However, since the neutrino-nucleon cross-section rises 
with increasing neutrino energy, this effect becomes non-negligible at neutrino energies 
above 10 TeV. 

The neutrino effective area Af5 is defined as the ratio of the rate of reconstructed 
events and the incident neutrino flux. It is calculated using Monte Carlo methods and 
incorporates the attenuation of the neutrino flux during the passage of the Earth, the 
neutrino-nucleon cross-section, the range of the generated muon and the reconstruction 
and selection efficiencies. This effective area is usually provided by the experimental 
collaborations. The energy dependence of the neutrino effective area comes mainly from 
the energy dependence of the cross-section (roughly oc E,) and the increase of the muon 
range. Notice that the muon effective area A^^, on the other hand, is defined as the ratio 
of the rate of reconstructed events and the incident muon flux. This area incorporates 
only the geometry of the detector and the detection efficiency. It is roughly equal to the 
geometrical area but might have a slight energy dependence. 

For the calculation of the spectrum of muon neutrino induced muons at the detector 
position we have to take into account the shift to lower energies due to the energy loss 
during muon propagation through matter |34] : 

oo RiJ.{Ev^,Efj) 



dE„ 



dn j dE,^, I dr e^'' 



dJ{E,,,e, 

dE,,.. 



{p — n) 



(3.10) 



where we neglected the attenuation term. In this expression the initial muon energy 
enters as an explicit function of the final muon energy as given by equation (13.81) . 
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Figure 4: Muon fluxes for tlie different decay channels of a dark matter candidate com- 
pared to the atmospheric background for upward through-going muons in standard rock. 
The flux is computed for a dark matter mass of 1 TeV (left) or lOTeV (right) and a 
lifetime of 10^^ s, for the neutrino spectra in Figures [2] and |3l In this case the muons lose 
energy on their way to the detector, smoothing out the spectral edges. 

Using equation flS.lOp we calculate the flux of through-going muons induced by neu- 
trinos from various dark matter decay channels and show the results in Figure H] for the 
case of a detector surrounded by standard rock. However, the result is also applicable 
for the case of detectors surrounded by water or ice since the dependence on the density 
cancels in equation fl3.10p and the muon energy loss parameters are roughly similar for 
the different materials ( cf. Table [1]) . Since there is no possibility to veto for the over- 
whelming background of atmospheric muons, only upgoing events and therefore a solid 
angle of 27i can be used for the analysis. We see that the deep inelastic scattering trans- 
forms the monochromatic neutrino lines into a continuous muon spectrum. In addition, 
the energy loss in the muon propagation smooths out all the spectra making the edge 
corresponding to half the dark matter particle mass less clear. Still the spectrum for a 
line signal remains steeper than the others at the endpoint. 



3.2.2 Contained Muons 



These events are similar to through-going muons but in this case the neutrino-nucleon 
interaction takes place inside the instrumented volume. If the muon track ends inside the 
detector the events are called contained. If the muon track leaves the detector one speaks 
of a partially contained event. The rate of muon neutrino induced (partially) contained 
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Figure 5: Spectra of contained muons for the different decay channels of a dark matter 
candidate compared to the atmospheric background. The event rate per km^ of detector 
volume (filled with ice) is computed for a dark matter mass of 1 TeV {left) or 10 TeV 
{right) and a lifetime of 10^^ s, for the neutrino spectra in Figures |2] and [31 The line 
signal is changed into a muon continuum due to deep inelastic scattering, but it retains 
a hard edge at half the decaying particle mass. Also the other spectra are softer than 
the original neutrino ones, with the continuum neutrino spectra producing a rise at low 
energy, unfortunately well below the background. 



track-like events per unit detector volume is given by 



- fan fdE. 



dE^dVdt J J ^ dE^^ 



— np + (p -)■ n) 



dE^ 



(3.11) 

where we also neglected the attenuation term. In this case also the hadronic cascade is 
contained in the detector volume and therefore, by measuring the energy of the muon as 
well as of the hadronic cascade, it is in principle possible to reconstruct the total energy 
of the incident muon neutrino. In this case, however, one has to rely on the detection 
also of the hadronic cascade which, as we will discuss later, seems to be challenging. 

The effective volume of the detector for contained events corresponds roughly to the 
geometrical volume (apart from boundary effects and reconstruction efficiency) and it 
is not enhanced by the muon range, which as we have seen, grows as E^. Therefore, 
the statistics for contained events is much lower than for through-going events at large 
energies. For instance in the case of Super-Kamiokande the event rate above roughly 
10 GeV is dominated by through-going muons. On the other hand, in the energy range 
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of interest for dark matter searches the muon range is of the order of one kilometer and 
therefore the expected rate of contained muons is comparable to the rate of through- 
going muons in detectors of cubic kilometer size. Thus, these contained events might be 
equally important for dark matter searches at the new generation of neutrino telescopes. 

In addition, for downgoing contained muon events there is the interesting possibility 
to reduce the background of atmospheric muon neutrinos by the detection of a coincident 
muon that was produced in the same parent meson decay [35]. This strategy could be used 
to increase the signal-to-background ratio for this channel, especially at large energies. 
However, we will not discuss this strategy quantitatively in this work. 

In Figure [5] we show the muon spectra for contained events calculated using equa- 
tion (13. lip for the case of a detector volume filled with ice. The result for a volume of 
water can easily be obtained rescaling the rate with the slightly different density. In this 
case there is no smoothing due to muon energy loss as in the case of through-going muons 
and the edges of the spectra are clearer, in particular for the case of a two-body decay. 

Here we only discussed the case where only the muon is measured since this is what 
can be done by the experiments at the moment. If the hadronic shower is also measured 
the combined reconstructed spectra would be as in Figures [2] and [3l This is similar to 
the case of electron and tau neutrinos that is discussed in the next section. However, 
as will be discussed there, that channel offers a better signal-to-background ratio and a 
better energy resolution and will therefore be of more interest once the showers can be 
measured and used for analyses. 

3.3 Electron and Tau Neutrinos 

The charged- current deep inelastic scattering of an electron neutrino off a nucleus pro- 
duces a hadronic shower and an electron that immediately causes an electromagnetic 
shower. The charged- current deep inelastic scattering of a tau neutrino off a nucleus 
produces a hadronic shower and a tau lepton. Due to the short lifetime of the tau lep- 
ton {cTr = 87.11 yum), at these energies it decays almost instantly and produces another 
shower at the interaction point. Thus, at energies below many TeV, detectors like IceCube 
cannot distinguish electron neutrino from tau neutrino events since both types produce 
similar showers in the detector [36]. In these cases, however, the whole neutrino energy 
is deposited in the detector and therefore it may be possible in principle to reconstruct 
better the initial neutrino spectrum. On the other hand, the analysis for cascade-like 
events is much more difficult than the analysis for muon tracks. No cascade events from 
atmospheric neutrinos have been identified yet and there are only first studies on this 
topic e.g. by the IceCube collaboration [3?]. For this reason there is no effective area for 
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Figure 6: Spectra of cascade-like events for the different decay channels of a dark matter 
candidate compared to the atmospheric background. The event rate per km^ of detector 
volume (filled with ice) is computed for a dark matter mass of 1 TeV (left) or lOTeV 
(right) and a lifetime of 10^^ s, for the neutrino spectra in Figures [2] and [3l 



this type of events available yet and therefore it is difficult to estimate realistically the 
sensitivity in shower events. 

Shower-like events are also characteristic of the neutrino-nucleon neutral-current 
interaction and for this reason probably only a combined analysis of neutral-current 
interactions for all neutrino flavours and charged-current interactions for tau and electron 
neutrinos will be feasible. In this case the total rate of neutrino- induced shower-like events 
is given by 



dN 



shower 



dVdt 



dQ 



E 



-e,T 
oo 



dJuiiE^^, 9, ' 



dE, 



[ac^E^^] Up + {p ^ n)] 



i^l -'-'shower 



+ E j ^^-^ 

-^shower 



dJui,{E, 



dE, 



ui ; -^shower ) 



dE. 



shower 



(3.12) 



We give in Figure [H] the signal and atmospheric background spectra calculated from 
equation fl3.12p for the case of a detector volume filled with ice. Note that in this case 
the muon neutrinos contribute only via neutral-current interactions which are weaker 
by a factor of about three compared to charged-current interactions (c/. equation ( 13. 5p ). 
Still, since the atmospheric muon neutrino flux is a factor of 20 larger than the electron 
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neutrino flux at TeV energies, the atmospheric muon neutrinos provide the dominant 
background. At the same time the signal is increased by roughly a factor of three. This is 
because, due to neutrino oscillations, the signal is roughly equal in all neutrino flavours 
and, therefore, the signal rate from the charged-current interactions of electron and tau 
neutrinos is the same as for the muon neutrinos. In addition, the combined neutral- 
current signal of all flavours contributes at the same level as the charged-current signal 
of one flavour. In summary, cascade-like events will offer a signal-to-background ratio 
that is roughly one order of magnitude larger than in the muon case and, therefore, 
they appear to be a very promising channel, if they are measured. We see also that in 
this case, assuming that the total shower energy can be reconstructed, the line-feature 
is preserved and clearly visible. 

4 Rates and Bounds 
4.1 Super-Kamiokande 

Super-Kamiokande is a 50 kt water Cherenkov detector. The fiducial mass is 22.5 kt and 
the muon effective area is 1200 m^ (with a slight zenith angle dependence due to the 
cylindrical shape of the detector). It is identical to the geometrical area since the re- 
construction and selection efficiencies are virtually 100 %. Super-Kamiokande has been 
looking for a neutrino signal, mostly from dark matter annihilation in the centre of the 
Sun, the centre of the Earth and in the Galactic centre. No excess has been found so far, 
and this can also be used to put a constraint on the decaying dark matter case. We com- 
pare the fiux of upward through-going muons from dark matter decay (integrated over 
energies above the threshold at 1.6 GeV) with the 90 % CL. fiux limit of excess neutrino- 
induced upward through-going muons provided by the Super-Kamiokande collaboration 
for the Galactic centre (the limits from the Sun and Earth fiux are weaker) [3H]. As 
discussed in section 12.21 the strongest bounds are obtained for the largest field of view. 
The exclusion region in the decaying dark matter parameter space derived from the limit 
in the 30° half-angle cone around the Galactic centre is given in Figure [71 

The bounds obtained here become stronger for larger masses although the neutrino 
fiux is proportional to I/ttidm for constant lifetime. This is due to the increasing neutrino- 
nucleon cross-section and the increasing muon range. The bounds are stronger for the 
two-body decay signal compared to the other cases since there the signal is concentrated 
at the end of the spectrum and benefits from the larger neutrino-nucleon cross-section 
and muon range. Note that these present bounds do not have sensitivity to the parameter 
region preferred by the PAMELA excess yet, which corresponds to a lifetime of the order 
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Figure 7: 90 % CL. exclusion region in the lifetime vs. mass plane for a decaying dark 
matter candidate from the non-observation of an excess in the Super-K data. The bound 
is stronger for a line signal, since there the spectrum is harder, resulting in a larger 
muon flux due to the increasing neutrino-nucleon cross-section and muon range. For the 
channels that contain Z° or bosons in the final state the exclusion range is cut at 
the threshold for their production. 

of 10^^ s and masses larger than 200 GeV. 

4.2 Rates and Bounds for Present and Future Experiments 

Assuming decaying dark matter with a lifetime of 10^^ s, we can now compute the 
expected signal rates for present and future experiments. These results can be easily 
generalised to arbitrary lifetimes, by recalling that the flux is proportional to I/tdm- 
We give the rates for some typical detectors of different sizes, i.e. Super-Kamiokande, 
ANTARES/AMANDA and IceCube. The results for Super-K can be easily scaled up to 
the Hyper-Kamiokande/UNO size by multiplying by a factor 10 or 20 (for a Hyper-K 
mass of 500 kt and Hyper-K/UNO mass of 1 Mt, respectively). The result for KM3NeT 
will be very similar to that expected for IceCube. 

We would like to stress here that Super-K is still taking data, and that the 
full ANTARES detector was completed in summer 2008 and is also operational. The 
AMANDA detector was decommissioned in summer 2009, but has since been substi- 
tuted by the partial IceCube detector, which already had 59 strings deployed in the 
ice in early 2009. The other experiments are still in the planning phase: KM3NeT is a 
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proposed cubic-kilometer sized underwater neutrino telescope in the Mediterranean Sea, 
which will probably have an effective volume comparable to IceCube, but will be able to 
look at the Galactic centre, while the proposed Hyper-Kamiokande and Underwater Neu- 
trino Observatory (UNO) are water Cherenkov detectors similar to Super-Kamiokande 
but of megaton scale. 

For the case of IceCube we also take into account the DeepCore subdetector, which 
is currently under construction. It is designed to lower the energy threshold of the ex- 
periment to roughly 30 GeV and to increase the sensitivity at low energies. This detector 
consists of six additional strings with less spacing between the digital optical modules 
compared to IceCube. One of these strings has already been deployed in early 2009. 
Commissioning of the full detector is planned for early 2010, while IceCube itself will 
only be completed in early 2011. The combination of IceCube and DeepCore can use the 
outer layers of IceCube as a veto to atmospheric muons and therefore has a An sensitivity 
for fully and partially contained events, but a considerably smaller effective volume. 

For the calculation of rates of upward through-going muon events we use the neu- 
trino effective areas for AMANDA, ANTARES and the IceCube 80 strings configuration 
from [39] and integrate over the muon spectrum. For the combined IceCube -|- DeepCore 
detector we amend the effective area in the low-energy range using the neutrino effective 
area given in [IQ]. In the case of Super-Kamiokande we calculate the rate using equa- 
tion (13. 9p with standard rock as surrounding material, a muon effective area of 1200 m^ 
and a threshold muon energy of 1.6 GeV. 

We see from Tables |2] and [3] that a sizable number of events is expected for a lifetime 
of 10^^ s, especially for experiments of cubic kilometer scale, such as to become significant 
above the atmospheric background even for a dark matter particle mass of 300 GeV. Of 
course for larger masses the significance becomes greater due to the increasing signal 
rate. Note that here we did not make use of any spectral information. In that case larger 
dark matter masses would also benefit from the falling background. 

Requiring the combined number of signal and background events not to exceed the 
background above the 90 % C.L. (in the Gaussian approximation this corresponds to 
a = S/\fB < 1.28), similar to the case of Super-K in Figure [71 we can then give in 
Figure |8] a forecast of the exclusion region which may be obtained from kilometer-cubed 
experiments using one year of data. The larger statistics of the future experiments will 
improve the Super-K bounds by more than an order of magnitude and explore the region 
of lifetimes above 10^^ s, for masses larger than 200 GeV. Note that for ten years of data 
the lifetime limit will become stronger approximately by a factor of three. For lower 
masses the bounds will remain weaker, but in that parameter region a very important role 
will be played by DeepCore, which will considerably improve the IceCube performance 
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Table 2: Number of upward through-going muon events per year from the atmospheric 
neutrino background and different dark matter decay channels at several neutrino ex- 
periments. The signals are given for a dark matter lifetime of 102^s and a dark matter 
mass of 300 GeV. 



between 30-100 GeV masses as can be seen in the right panel of Figure [8], and also by 
the megaton water detectors which are expected to strengthen the Super-K bounds by 
an order of magnitude down to masses of a few GeV. This low-mass region does not 
provide an explanation of the PAMELA excess and is plagued by a stronger atmospheric 
background, but it is still remarkable that even there lifetimes larger than 102^ s will be 
probed in future experiments. 



4.3 Energy Resolution and Reconstructed Spectra 

Once a signal has been detected, the question arises if it will also be possible to re- 
construct the neutrino spectra and extract some information on the dark matter decay 
channel. For this purpose one important factor is the energy resolution of the neutrino 
detectors. We will take here for reference the IceCube detector, for which the energy res- 
olution is logio(-Emax/-E'min) = 0.3-0.4 for track-like events and logio(-E'max/-E'mm) = 0.18 
for cascade-like events [H]. 

We show in Figures [914T41 the histograms for the signal and the atmospheric back- 
ground using an energy resolution of 0.3 in log^^g ^ three bins per decade for upward 
through-going and contained muons, and an energy resolution of 0.18 in log^o-E' and 
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Table 3: Number of upward through-going muon events per year from the atmospheric 
neutrino background and different dark matter decay channels at several neutrino exper- 
iments. The signals are given for a dark matter lifetime of 102^ s and dark matter masses 
of ITeV (top) and lOTeV (bottom). 
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Figure 8: 90 % CL. exclusion prospects in the lifetime vs. mass plane for a decaying dark 
matter candidate from the non-observation of a significant excess in the total rate of 
neutrino induced upward through-going muons observed at IceCube ileji) or IceCube -|- 
DeepCore (right) in one year. Clearly visible is the enhanced sensitivity of the DeepCore 
extension in the low-mass region. 

five bins per decade for shower events. These figures can be compared to Figures |4H6] 
which show the spectra unbinned and without finite energy resolution. Also shown is the 
significance of the signal over the background in different bins for a lifetime of 10^^ s for 
the different channels using one year of data with an effective area of 1 km^ for upward 
through-going muons and an effective volume of 1 km^ for contained muons and cascades. 
All plots are available for both, scalar and fermionic, dark matter candidates. 

We see that for nearly any of the spectra, the signal will appear with a large statistical 
significance in more than one single bin and it will be clear that the neutrino signal is 
not following a power law like the atmospheric one. Thus, it is clear that using spectral 
information it will be possible to set much stricter limits on the decaying dark matter 
parameter space than shown in Figure [H In order to give an idea of the sensitivities 
that can be obtained using spectral information, we show in Table H] the values of the 
dark matter lifetime for several decay channels that correspond to a 5 a signal in the 
most significant energy bin after one year of observation for an idealised detector with an 
effective muon area of 1 km^ and an effective volume of 1 km'^ for contained muons and 
shower events. We see there that the limits from through-going and contained muons 
are better but not far from those shown in Figure [HI while the shower events in principle 
allow to reach even one order of magnitude larger lifetimes. Using not only the dominant 



24 



DM 

DM (T-u) 
DM ^ ZZ (WW) 
atmospheric 




DM ^ vv 
DM (tt) 
DM -> ZZ (WW) 



through-going muons 
moM = 1 TeV 



CO 

II 



'•DM 



10° 



10'' 



100 1000 10000 



(GeV) 



10" 

35 
30 
25 
20 



DM vv 
DM (tu) 
DM ^ ZZ (WW) 
atmospheric 




through-going muons 



mDM = 10TeV, Tdm = 10^^s !- 



10 100 1000 10000 

(GeV) 



C/3 

II 

t) 15 
10 
5 



DM vv 
DM ^ (1(1 (tt) 
DM ^ ZZ (WW) 



through-going muons 
mDM = 10TeV 

"■•DM = 10^^ S 



10 



1 00 1 000 
E^ (GeV) 



10000 




1 00 1 000 
E^ (GeV) 



10000 



Figure 9: Top: Muon fluxes for tlie different decay cliannels of a scalar dark matter can- 
didate compared to the atmospheric background for upward through-going muons. The 
flux is computed for a dark matter mass of 1 TeV (left) or 10 TeV (right) and a lifetime 
of 10^^ s using an energy resolution of 0.3 in log^Qi? and three bins per decade. Bottom: 
Statistical significance of the signal of through-going muons shown above calculated for 
every single energy bin using one year of data with an effective area of 1 km^. 
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Figure 10: Top: Muon rates per km^ of detector volume for the different decay channels 
of a scalar dark matter candidate compared to the atmospheric background for contained 
muons. The flux is computed for a dark matter mass of 1 TeV (left) or 10 TeV (right) and 
a lifetime of 10^^ s using an energy resolution of 0.3 in log^Q E and three bins per decade. 
Bottom: Statistical significance of the contained muon signal shown above calculated for 
every single energy bin using one year of data with an effective volume of 1 km^. 
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Figure 11: Top: Shower rates per km'^ of detector volume for the different decay chan- 
nels of a scalar dark matter candidate compared to the atmospheric background for 
electromagnetic and hadronic showers. The flux is computed for a dark matter mass 
of 1 TeV [left) or 10 TeV (right) and a lifetime of 10^^ s using an energy resolution of 
0.18 in log^o ^^"i fi'^^ bins per decade. Bottom: Statistical significance of the shower 
signal shown above calculated for every single energy bin using one year of data with an 
effective volume of Ikm^. 
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Figure 12: Top: Muon fluxes for tlie different decay channels of a fermionic dark matter 
candidate compared to the atmospheric background for upward through-going muons. 
The flux is computed for a dark matter mass of 1 TeV (left) or 10 TeV (right) and a 
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Bottom: Statistical significance of the signal of through-going muons shown above calcu- 
lated for every single energy bin using one year of data with an effective area of 1 km^. 



28 



10' 



10° 



;io" 



DM 



DM^Zv 
eev ((ihv/ttv) 
DM ^ We 
DM ^ (Wx) 
atmospheric 




100 1000 
(GeV) 



10000 



DM ^Zv 

DM ^ eev (hhv/ttv 
DM We 
DM Wn (Wx) 



contained muons 



m 



DM 



= 1 TeV 
lO^^s 



10 



100 1000 
(GeV) 



10000 



10' 



10° 



DM 



DM^Zv 
—> eev (hhv/ttv) 
DM ^ We 
DM ^ Wn (Wx) 
atmosplieric 




100 1000 
E^ (GeV) 



10000 





18 




16 




14 




12 


m 






10 


CO 




II 

D 


8 




6 




4 




2 








DM ^Zv 

DM eev (|j,(iv/xxv) 
DM ^ We 
DM ^ Wn (Wx) 



contained muons 
mDM = 10TeV 



10^6 s 



10 



1 00 1 000 
E^ (GeV) 



10000 



Figure 13: Top: Muon rates per km'^ of detector volume for the different decay chan- 
nels of a fermionic dark matter candidate compared to the atmospheric background for 
contained muons. The flux is computed for a dark matter mass of 1 TeV (left) or 10 TeV 
{right) and a lifetime of 10^^ s using an energy resolution of 0.3 in log^g E and three bins 
per decade. Bottom: Statistical significance of the contained muon signal shown above 
calculated for every single energy bin using one year of data with an effective volume of 
1 km^. 
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Figure 14: Top: Shower rates per km'^ of detector volume for the different decay chan- 
nels of a fermionic dark matter candidate compared to the atmospheric background for 
electromagnetic and hadronic showers. The flux is computed for a dark matter mass 
of 1 TeV (left) or 10 TeV (right) and a lifetime of 10^^ s using an energy resolution of 
0.18 in log^o ^^"i fi'^^ bins per decade. Bottom: Statistical significance of the shower 
signal shown above calculated for every single energy bin using one year of data with an 
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Table 4: Dark matter lifetimes corresponding to a cr = 5 significance in the most signif- 
icant energy bin after one year of observation in an idealised detector with an effective 
muon area of 1 km^ and an effective volume for contained muons and showers of 1 km^. 
The numbers are given for dark matter masses of 1 TeV (top) and lOTeV (bottom). 
Notice that the sensitivity obtained with through-going and contained muons is similar. 
At larger masses the bound from through-going muons is stronger since the statistics 
increases due to the longer muon range at higher energies. However, neglecting recon- 
struction efficiencies the strongest constraint is obtained from shower events since that 
channel offers the best signal-to-background ratio (see discussion in section l3T3ll . 

energy bin from Figures [9]-[T4l but a combination of several energy bins optimised for 
each individual decay channel it will be possible to set even stronger constraints on the 
dark matter lifetime. So a signal in the region preferred by PAMELA should be in the 
detectable range. 

On the other hand, discriminating between the spectra for the different channels will 
not be so straightforward, especially if the mass of the decaying particle is unknown. After 
convolution with the energy resolution, the two-body, three-body or continuum spectra 
appear quite similar, especially within their statistical error, but their significance peaks 
at slightly different values for the same dark matter mass. Note, though, that the signal 
from a neutrino line remains steeper than the other ones at the edge and it may be 
possible to distinguish it with sufficient statistics. In this respect the more promising 
strategy is probably exploiting the better energy resolution of the shower events, if they 
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can be detected. In general a comparison between the different types of events, through- 
going, contained and cascade-hke, will make disentangling the shape of the spectra easier. 
Moreover, if the dark matter mass is measured via another channel, like gamma-rays, it 
may be possible to exploit this information in the neutrino fit and compare the position 
of the neutrino signal "peak" in the data with the expectation. This should help in 
disentangling at least a continuum spectrum from the two- and perhaps also three- 
body decay cases. For this specific strategy probably one of the most promising dark 
matter candidates would be a fermion, which may decay into a leptonic three-body final 
state and subdominantly into 7Z/, as e.g. gravitino dark matter with trilinear i?-parity 
breaking |12]. Then, if the 7-channel is suppressed by 2-3 orders of magnitude compared 
to a leptonic three-body decay, the gamma-ray and neutrino experiments will actually 
be exploring the same range of lifetimes. In this case an observation of a gamma line by 
Fermi will provide the dark matter mass measurement and the neutrino signal with a 
much shorter lifetime will point at a three-body or Z^v dominant decay. In case both 
signals in gamma-rays and neutrinos are measured, it may be possible to disentangle 
also between a scalar and a fermionic dark matter candidate, which seems to be very 
difficult from neutrino measurements alone, since the two types of particles produce very 
similar spectra within the energy resolution of the detectors, as can be seen comparing 
Figures El [JO] and [H] with Figures [121 US] and [H 



5 Conclusions 

We have studied in this paper the possible neutrino signals from decaying dark matter, 
considering different decay channels and spectra, both for a scalar and a fermionic can- 
didate. We have concentrated here on the region of parameter space that is preferred 
in order to explain the PAMELA positron excess and shown that in this case a signal 
may soon be visible at neutrino observatories. The non-observation of such a signal will 
put rather strong constraints on the leptophilic decaying dark matter explanation of the 
excess, except for the decay into e+e~ where no signal in neutrinos is expected. In this 
sense neutrino observations are complementary to other astrophysical constraints com- 
ing from radio frequencies and inverse Compton emission, which are more sensitive to 
the electron channel. A neutrino signal will allow to disentangle between decaying and 
annihilating dark matter, by comparing the signal towards and away from the Galactic 
centre [19], and also between dark matter and astrophysical sources. 

More difficult is the identification of the dark matter decay channels, since all neu- 
trino spectra finally result in a broad peak in the muon spectrum. However, the analysis 
of cascade-like events, which contain in principle all the neutrino energy and have the ad- 
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vantage of a better energy resolution, may improve the situation and allow with sufficient 
statistics to disentangle at least a line-like feature. Also, the neutrino signal alone can- 
not distinguish between scalar and fermionic dark matter candidates since the resulting 
spectra are very similar in the two cases. 

On the other hand, for some of the decay channels discussed here, like the ones with a 
Z^/W^ gauge boson in the final state, corresponding signals are also expected in gamma 
rays and antiprotons and may provide additional information on the model parameters 
and a cross-check of the decay channel. Even for the pure leptophilic channels, gamma 
rays from final state radiation or from subdominant decays may play an important role 
in discriminating between models due to the better sensitivity in the gamma-ray channel. 
In general, neutrino observations offer complementary information and can be used to 
test also models where the gamma-ray signal is strongly suppressed compared to the 
leptonic one. The next generation of neutrino experiments can therefore be expected to 
yield some very interesting results. 

Note Added 

During the completion of this work the preprint [l3] appeared, presenting a related 
analysis. 
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